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Macroscopic quantum tunneling (MQT) for a single fluxon moving along a long 
Josephson junction is studied theoretically. To introduce a fluxon-pinning force, we 
consider inhomogeneities made by modifying thickness of an insulating layer locally. 
Two different situations are studied: one is the quantum tunneling from a metastable 
state caused by a single inhomogeneity, and the other is the quantum tunneling in 
a two-state system made by two inhomogeneities. In the quantum tunneling from a 
metastable state, the decay rate is estimated within the WKB approximation. Dis- 
sipation effects on a fluxon dynamics are taken into account by the Caldeira-Leggett 
theory. We propose a device to observe quantum tunneling of a fluxon experimen- 
tally. Required experimental resolutions to observe MQT of a fluxon seem attainable 
within the presently available micro-fabrication technique. For the two-state system, 
we study quantum resonance between two stable states, i.e., macroscopic quantum 
coherence (MQC). From the estimate for dissipation coefficients due to quasiparticle 
tunneling, the observation of MQC appears to be possible within the Caldeira-Leggett 
theory. 
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§1. Introduction 

It is commonly recognized that the sine-Gordon equation plays an outstanding role in 
many physical problems. One of the most important applications is a long Josephson junc- 
tion (LJJ). When the junction width is taken large enough in one direction (defined as 
x-direction), a phase difference of superconductors across the junction may have spatial 
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dependence in the x-direction. It is believed that dynamics of the phase difference is well 
described by a classical equationS^ 

(ptt - (pxx + sin0 + a(f)t - Pcpxxt + / = 0. (1.1) 

Here x and t are measured in units of the Josephson penetration length Aj and of the inverse 
Josephson plasma frequency u!~^, respectively. The dissipation coefficient a is related to 
quasiparticle tunneling through the oxide barrier, and (3 is related to the normal current of 
quasiparticle parallel to the junction. The external current / is assumed spatially uniform. 

It is known that solitons in the form of fluxons propagate along the junction following the 
classical equation ( p..l|) . Experimentally, fluxons in LJJs were first observed indirectly by 
zero-field steps on the current voltage (I-V) characteristics of the junction.!^ The zero-field 
steps are well explained by fluxon propagation governed by eq. ( |1.1| ) with repeated reflections 
at the open ends of the junction.!^ Since then, development of experimental techniques has 
made it possible to directly observe profiles of separate fluxons and a space-time pattern of 
their interaction.S & 

These experiments are basically explained by 'classical' theories, and there is no clear 
experimental indication so far which needs 'quantum' theories to explain its results. One 
might argue that it is natural to expect an essentially classical motion for a single fluxon 
of the size of micrometer. However, we show in this paper that quantum effects can indeed 
show up in this extended object. Recently, Hermon et al. discussed the quantum dynamics 
of a single fluxon in a long circular Josephson junction.!^ In a subsequent paper,& they 
discussed a possible fluxon interference experiment. However in those papers, decoherence 
effects due to couplings to the environment, i.e., dissipation has been neglected. From the 
study in dissipative two-level systems,^'' it is expected that the quantum effects are strongly 
suppressed by the dissipation so that in real experiments, the observation of the quantum 
effects proposed by Hermon et al. may be more difficult than in the ideal case without the 
dissipation. 

In this paper, we propose other experiments to observe quantum effects of a single fluxon. 
The effect which we deal with here is quantum tunneling of a fluxon. Quantum tunneling 
phenomenon is one of the most typical ones which cannot be explained by classical theories. 
Because a fluxon is a macroscopic object with a length scale of micrometer, the quantum 
tunneling of a fluxon can be recognized as the macroscopic quantum tunneling (MQT). The 
macroscopic tunneling phenomena have already been studied in some other systems. For 
example, much attention has been paid on MQT of the phase of current-biased junctions and 
the dissipation on it. MQT in this system has been observed experimentallyJi'0'' and it has 
been claimed that the tunneling rate agrees with the value predicted by the Caldeira-Leggett 
theory0 within a phenomenological treatment of the dissipation. 
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Fig. 1. Schematic drawings of (a) a microresister (jc > j'c) and (b) a microshort {jc < j'^). The normalized 
strength of pinning is given by e = b[jc — jc)/jcAj in both cases. Here b is the length of the area where 
the critical ern'rent is modified. 



To consider the quantum tunnehng of a fluxon, we introduce structural inhomogeneities 
which capture a fluxon at a certain point in a Josephson junction line. Such inhomogeneities 
have been studied by McLaughlin and Scott in a pioneering paper.@ In that paper, they 
have considered interaction between a fluxon and a microresister (i.e., a narrow region where 
critical current density jc is reduced to j^) or a microshort (i.e., an narrow region with a 
enhanced critical current density). Schematic drawings of these inhomogeneities are shown 
in Fig. 1. McLaughlin and Scott have proposed a model of LJJs with a local inhomogeneity 
in the form of the equation of motion 

0ti - (pxx + sin + a(j)t - P4>xxt + f -e6 (x) sin = 0. (1.2) 

The last term in the left hand side represents a local change of the critical current density 
at X = 0. The normalized strength of the inhomogeneity, e is deflned hj e = {jc — jc)b/jc^3- 
The cases e > and e < correspond to a microresister and a microshort, respectively. As 
we show later, when e > (i.e. a microresister), the fluxon is attracted to the microresister 
and captured there if the fluxon does not have enough kinetic energy. In other words, the 
microresister plays a role of a pinning potential for the fluxon. We only deal with the case 
e > 0, i.e., microresisters. Throughout this paper, we assume 6 <^ Aj (i.e. £ <^ 1) so that the 
soliton size is larger than b to justify the description by 6{x) in ( |L2|) for the microresister. 

In this paper, we study the quantum tunneling of a fluxon in two situations. First, we 
consider the situation that there exists only one microresister in the junction. When the 
external current / is applied on the junction, the fluxon-pinning state becomes metastable. 
We study quantum tunneling of a single fluxon from the metastable state. The decay rate 
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is estimated within the framework of the WKB approximation by taking proper values of 
experimental parameters. We also propose experimental design to observe quantum tun- 
neling of a fluxon, and study conditions to allow the observation in this device. Next, we 
study the quantum tunneling between two stable states made by two microresisters. In 
this situation, we pursue possibility to observe resonance between two localized levels, i.e., 
macroscopic quantum coherence. In both cases, we consider dissipation effects using the 
Caldeira-Leggett theory. 

The paper is organized as follows. In § ^ we introduce a model Hamiltonian for fluxons in 
the analogy to the Caldeira-Leggett theory, and derive the effective action which describes 
dynamics of a fluxon. In § |^, we summarize controllable experimental parameters, and 
set up them to the values which are accessible in practical experiments. Some quantities 
characteristic of LJJs, are also estimated. In § ^, we consider quantum tunneling of a single 
fluxon from a metastable state, and in § |^, we study quantum tunneling of a single fluxon 
in a two-state system made by two microresisters. Concluding remarks are given in § ^. 

§2. Formulation 

2. 1 Classical equation of motion for a single fluxon 

In this subsection, we derive a classical equation of motion for a fluxon. For this purpose, 
we analyze the equation ( [1.2D with the assumption that parameters, a, (3, f and e are all 
small. In this case, eq. ( |1.2|) can be considered as a perturbed sine-Gordon equation. The 
soliton solution of the unperturbed sine-Gordon equation takes the form of a kink. 



which corresponds to a fluxon in LJJs. Here q{t) = ut is the center coordinate of the kink, 
and u is the kink velocity < 1). The velocity is normalized by the light velocity, c = AjcUp. 
In this paper, we only discuss the nonrelativistic limit \u\ ^ 1. 

Based on the kink solution the classical perturbation theory can be applied to 

eq. ( |L2| ). McLaughlin and Scott have shown that the perturbations only affect dynamics of 
the center coordinate q{t) and do not change the form of the kink within the framework of 
the lowest approximation. 0^ The equation of motion for q(t) is obtained as 



in the nonrelativistic limit. Here m is the classical soliton mass and identically equals to 8. 
Note that the damping strength working on the dynamics of a fluxon is related to a and (3 






(2.2) 



which are the coefficients of the dissipative terms in the field equation ( |1.2| ) for (j){x,t). The 
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effective potential V{q) for a fluxon is given by 

V{q) = -2nfq - (2.3) 
cosri q 

The first term is a driving force due to tlie external current, and the second term is a pinning 
potential caused by the microresister located at the origin. 

2.2 The Caldeira-Leggett theory 

As shown in eq. ( |2.2| ), the classical equation of motion for a fluxon includes a damping 
term due to the fact that the phase difference is a macroscopic variable. To deal with 
dissipation effects on MQT of a fluxon phenomenologically, the Caldeira-Leggett mo IS 
introduced in this subsection. We provide harmonic oscillators coupled to the macroscopic 
variable q, and write the Hamiltonian as 
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:2.4) 



Here p = mq is a momentum variable conjugate to g, whereas qj and pj = rrijqj are a 
coordinate and a momentum of harmonic oscillators, respectively. The Caldeira-Leggett 
model has succeeded in explaining experiments of junctions phenomenolog ically.© Here, we 
extend this treatment to LJJs. The reservoir parameters, rrij, uj and cj are characterized 
by a spectral function 

^M = ^E^^(^-^.)- (2-5) 

2 Y ^j^j 



We choose the spectral function as 



J{u) =mla + ^juj. (2.6) 



Then, the classical equation ( p.2|) is derived from the Hamiltonian (|2.4| ) by eliminating the 
reservoir degrees of freedom.^ 

From this model, we obtain the partition function Z after integrating out the reservoir 
degrees of freedom as 

Z = Znf Vq{r) exp (^-^S,f, [ q{r) ]^ . (2.7) 

Here is the partition function of the unperturbed reservoir, and g"^ is the normalized 
Planck constant defined by g^ = Uuo^/Eq where Eq is the energy scale in the LJJ. The 
effective action Scs is calculated 

S.Mr)] = fj^ drr-me + V{q)\ + \ fj"^ dr dr'Kir - r') Ur) - g(r')] , (2.8) 
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where T is a temperature normalized by hup/k^. Here, is Boltzmann's constant. The 
kernel K{t) is given by 



cu(l/2r- \t\ 



cosh 



(2.9) 



(2.10) 



sinh(tu/2T) 

In the later sections, we make estimate of the tunneling rate based on the effective action 



The derivation of the effective action (p.8|) is, however, quite intuitive. In the above 
derivation, infinite degrees of freedom of the field (p is first reduced to only one degree of 
freedom, i.e., the center coordinate q in the equation of motion ( p.2|) . Then, quantization 
for q is performed. If the quantum fluctuation around a fluxon affects quantum tunneling, 
the effective action (p.8|) may take a different form, because the fluctuation of the field is 



eliminated in the derivation of ( |2.2| ) before the quantization. In the next subsection, we 
check the validity of the effective action ( |2.8|) by applying the Caldeira-Leggett formalism 
directly to the equation of motion for the field variable, eq. ( |1.2| ). The readers who are not 
interested in details of theoretical derivation may skip the next subsection. 

2.3 Another derivation of the effective action 

In this subsection, the effective action (|2.8|) is derived in another way. The plan is as 
follows. First, quantization for the field variable is performed. The dissipation effects are 
taken into account by constructing a model Hamiltonian on the analogy of the Caldeira- 
Leggett theory. Then, the semiclassical theory is applied to the model with the path integral 
method. Through the perturbation expansion for the parameters, a, f3, f and e, we show 
that the effective action (p.8|) is reproduced within the lowest order approximation. 



First, we derive the Hamiltonian which yields classical field equation ( |1.2| ). In the absence 
of dissipation, the Hamiltonian is easily obtained as 



dx 



^0x + ^0? + (l-cos. 



+ /0-£<5(x)(l 



COS( 



(2.11) 



In the dissipative case, the dissipative terms, —acpt and Pcpxxt are obtained by introducing 
two kinds of harmonic oscillators coupled linearly via the field variables 0(a;, t) and (j)x{x,t) 
as follows: 



jdxY: 



2m i 



+ 



qj{x) 



mjUJj 



(2.12) 
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2m'j 2 



q'Ax) 
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/ /2 



(2.13) 



These oscillators constitute a proper heat bath causing the dissipation. 
The reservoir is characterized by the spectral functions Ja{^) and J/3{uj) defined by 

M^) = - y]—^S{^-^j), (2.14) 



2 Y "^i^i 



J2 



2 Y "^i^i 



(2.15) 



In order to produce the dissipative terms, —a(j)t and l3(j)xxt, we choose the spectral functions 
as 



Ja{uj) 



(3uj. 



The total Hamiltonian is given by 



H = Hs + Ha + Hg. 



(2.16) 
(2.17) 

(2.18) 



From this Hamiltonian, the perturbed sine-Gordon equation ( |1.2| ) is derived in the classical 
limit. The details of the derivation is given in Appendix [^. The dissipative sine-Gordon 
system with the Hamiltonian Ha has been used in the study of charge density wave0^ and 
long Josephson junctions.lii' It should be noted that the dissipation term /3(pxxt in ( P^ ) is 
expressed in a simple way by the coupling to the derivative 0a;- 

After integrating out the the reservoir degrees of freedom, we obtain the partition function 
Z as 



Z = ^P0exp Lls,s[^ix,T)] 



(2.19) 



Here Zj^'' and Z^' are partition functions of the unperturbed harmonic oscillations, and 
is the normalized Plank constant defined in § |2.2| . The effective action Ses is given by 



eflf 



5*8 + + 5*, 



/3' 



where 



Ss 



dr / dx 
1 /-VT 



+ ^0? + (1 -cos( 



- f4> — e5{x) (1 — cos 0)^ 
dr dr' j dxKair - r') ^0(x, r) - 0(x, r') j , 

Sp = \j^ dr^ dr' j dxKp{T - t'){^x{.x,t) - (j)x{x,T')^ . 



2 Jo 



(2.20) 

(2.21) 
(2.22) 

(2.23) 
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The kernels, Kair) and Kj^It) are given by 



1 r°° 

Ka{r) = - / aLuD^iT)du;, (2.24) 



TT Jo 
1 /""^ 



/ /3cuD^(r)da;, (2.25) 

JO 



TT Jo 

where D^{t) is defined by (2.10). 

Here we derive an effective theory of a single fluxon by approximating the quantization 
of the field theory defined by (2.21)-(2.23). For this purpose, we consider the semiclassical 
theory for the sine-Gordon equation. The semiclassical quantization of the field theory 
has been studied extensively in the context of high-energy physics in IQTO's.S'lll' In the 
semiclassical method, the classical soliton solution is regarded as the ground state of a Fock 
space called the one-soliton sector. This sector is completely separated from the sector 
containing no soliton, because of the topological stability of the soliton. The states of the 
one-soliton sector are constructed by a perturbative expansion in g"^. This expansion is valid 
as long as is small. As shown in § g"^ is estimated as g^ < 10~^ in proper choices 
of experimental parameters. Hence, the semiclassical approach can be applied. In this 
subsection, we only explain the outline of this approach. The details of the semiclassical 
calculation is given in Appendix [^. 

First we consider the unperturbed sine-Gordon theory. The partition function is given by 

Zo = fV(f){x,T)exp |^-l^o[0(x,r)]j . (2.26) 
Here 5*0 is the Euclidean action, 

So = jdr jdx {^4>^ + ^cg + (1 - COS0)) , (2.27) 

where (j) = tt- This field theory possesses a nontrivial stationary solution 
at 



(x — g) = 4 arctan 



exp(x — q) 



(2.28) 



where q is the center coordinate of the soliton. This new variable q is called a collective 
coordinate in the field theory, and regarded as a dynamical variable. In the semiclassical 
approach, we only consider the paths around the stationary solution because they predom- 
inantly contribute to the partition function (|2.26|) . The deviation around the stationary 
solution is denoted by rj as 

0(s, r) = 0o(s - g(r)) + gr,{x - q{r), r). (2.29) 

In terms of LJJs, the collective coordinate g(r) describes dynamics of a single fiuxon, and the 
rest degrees of freedom, r]{x, r) represent the fiuctuation around the fiuxon. The perturbation 



8 



expansion of the action ( p. 27 ) in terms of rj is nothing but that in terms of g"^. In the lowest 
approximation, we obtain 'plasmon' excitations by the quantization of the quadratic parts 
for 1] in the action. The plasmon excitations have a energy gap of hup. Since the energy 
gap is estimated as several Kelvin, the plasmon excitation can be neglected at low enough 
temperatures of the order of mK. The higher order terms of rj in the action only contribute 
to the partition function as the higher order of g^. The semiclassical expansion of g^ needs 
careful treatments especially in the treatment of the collective coordinate (See Appendix 
We only write the result: 

Zo^ f Vq{r) exp y dr (^m + ^ . (2.30) 

From this result, we find that the fluxon behaves like a free particle with the mass m = 8. 
It is a natural result because there is no force which modifies the velocity of the soliton. 



Next we consider the perturbed sine-Gordon equation (|1.1| ). The perturbation part of the 
action is given by 

5'ext[0(x,r)] = J^'^^dT Jdxf(P, (2.31) 

rl/T r 

Spi4(j){x,r)] = J dr Jdxe6{x){l-cos(f)), (2.32) 

in addition to (2.22) and (2.23). Here we assume that the perturbation parameters, a, (3, 
/, e are all small. Then, we take the perturbation expansion in terms of these small terms 
in addition to the expansion in g^. From straightforward calculation,^* it is shown that the 
lowest order contribution of the perturbation expansion is obtained only by substituting the 
soliton solution ( p.28| ) to the perturbative actions. In other words, the perturbations do not 
modify the waveform of the soliton in the lowest order. 
First, we consider the external current term (2.31). Substitution of (|2.28|) gives 

SLM^)] = '^oxt[0o(a; - g(r))] - 5'ext[0o(a;)] 

r-l/T 

dr(-27r/g). (2.33) 



Here the prime implies the effective action for the variable q{T). We subtracted the constant 
Sext[(po{x)], because we have chosen the origin of the potential energy at g = 0. 
The pinning potential is obtained by substituting (|2.28|) to (2.32), 

S'pinW)] = Spin [00 (a; - q{r))] 
cosh^ q 

Finally, we consider the dissipation described by (2.23) and (2.22). The effective action is 
obtained as 

S'Mr)] = S4Mx - qir))] 



d^TTTTT-- (2-34) 
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= ^/^'''^dr^"dr'y'da;ir„(r-r')(^0o(^-g(T))-0o(a;-g(r'))j , (2.35) 
S'Mr)] - S^iMx - q{r))] 

= ^/^'^^dT/^'dT'/dxi^^(r-rO[0[)(a;-g(r))-0;(a;-g(r'))) , (2.36) 



where 4>q = Since the integration over x is difficult to carry analytically, we assume 

that the paths of the collective coordinate g(r) satisfy 

|g(r)-5(r')|«l (2.37) 

for all T, t' . This condition is well satisfied in the case we consider in this paper. Using the 
approximations 

- q{T)) - - Qir')) ^ -<Poix - q{r)) (q{r) - g(r') j , (2.38) 
0o(^ - g(r)) - <P',{x - q{r')) ^ -^[((a; - q{r)) (^g(r) - g(r')) , (2.39) 
and the identities 

j dx0o = 8 = m and j dx0o^ = y, (2.40) 
we obtain from (2.35) and (2.36) 

S'Aq{r)] = \ l'^^ dr dr'mK^ir - r') [q^r) - g(r')) (2.41) 
S'M'')] = \ - U^^ - ^(^'))'- (2.42) 



We summarize the total effective action of a fluxon. The partition function is given by 

Z = fvq{r) exp (-lsMr)]\ (2.43) 



where 



S'Aqir)] = / dr (^me + V{q)) + ^ / dr / dr'X(r - r') (^?(r) - g'(r)) ' (2-44) 

V{q) = -27: fq - (2.45) 
cosh [q) 

The kernel K{t) is given by 

K(r) ^^Jm(a + (^^ D^doo, (2.46) 
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Fig. 2. A layout of a LJJ. The junction length L is taken as large enough compared with the junction 
width W. 



where D^^{t) is defined in eq. (2.10). The effective action obtained here agrees with one 



given in eq. (|2.8|) , and it has been shown that more intuitive derivation in § |2.2| is justified. 
One might feel that this result is trivial. It is, however, nontrivial because the fiuxon has a 
finite size in space. In fact, if the condition ( |2.37| ) does not hold, the action has a different 
form from the simple action (2.44). 

§3. Experimental Parameters 

In this section, we summarize controllable experimental parameters relevant to the quan- 
tum tunneling of a fiuxon. A layout of a LJJ is shown in Fig. 2, and experimentally param- 
eters realizable within the available technique are given in Table I. The values given in the 
table are indeed typical ones in actual experiment s.^-* The width of the Josephson junction 
W is now a controllable variable and scaled by micrometer. 

First, we associate the experimental parameters with the parameters in the normalized 
sine-Gordon equation ( |1.1|) . In the derivation of (|1 . 1| ) , the normalized length, time and 
energy are introduced. In the original unit, the equation of motion without inhomogeneities 
is written as 



0tt - Aj0^^ + sin0 + au^ ^(f)t - jSu^ ^\](t)xxt + 3 / 3c = 0. (3.1) 



Here, j is the external current density applied to the junction. 

The length scale in the LJJ is given by the Josephson penetration length Aj, which may 
be varied in the range from 20 [yum] to 200 [/im] by controlling the critical current density 
jc-^ In the present estimate, Aj is obtained as 

Aj = ( ^ ) ~ 27 [/iml. (3.2) 
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Table I. Experimental parameters used for an estimate. 



ParaiTK^tcTS 



Estimat(xl value 



d : Thickness of the junction 2 [nm] 

Sr : Relative dielectric constant 8 

jc : Critical current density 3 x 10^ [A/m^] 

W : Width of the junction W [/um] (variable) 

Xl ■ London penetration length 60 [nm] 



Here $o = h/2e is the unit flux and /xo(= 47r x 10~^ [Hm~^]) is the permeabihty of the 
vacuum. The size of a fluxon is characterized by Aj and the junction length L must be taken 
larger than Aj so that at least a single fluxon exists in the junction. 

The time is scaled by where cUp is the plasma frequency, which is estimated from the 
value given in Table / as 

where So{= 8.85 x 10~^^ [Fm~^]) is the dielectric constant of the vacuum. Note that Up is 
independent of W. From this estimate, the plasmon excitation gap is estimated as 

;ia;p~ 3.9 [K]. (3.4) 

Well below this temperature, the plasmon excitations can be neglected. 
The energy is measured by a unit energy Eq, 

Eo = ^jcWX, ~ 1.93 X lO^W [K]. (3.5) 

Prom this, the energy of a single fluxon is obtained, 

8Eo ~ 1.55 X lO^W [K] (3.6) 

where the factor 8 comes from the dimensionless mass of a fluxon. The energy of a single 
fluxon is so large that nucleation of soliton-antisoliton pairs can be neglected at sufficiently 
low temperatures. 
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The normalized Planck constant, g is given as 



1/2 



Note that is independent of the critical current jc- In the present estimate, the value of 
g^ becomes 

, 2.0 X 10-3 , , 

/ - (3-8) 

When W ~ l[yuni], we get small value, g^ ~ 0.002. Hence, the semiclassical approach 
explained in § ^]3| is applicable. 
The mass of a fluxon is estimated as 

mf = — ~ 1.25 X IQ'^W me, (3.9) 

where me is the electron mass and c is defined by c = Ajcjp. Since we find that the fluxon 
mass is remarkably small (~ lO-'^me if we take as 1 [/im]), we can expect substantial 
quantum effects of fluxons, in spite of the large fluxon size and the small normalized Planck 
constant. Moreover, it should be noted that mf is proportional to the junction width W . 
Hence, the fluxon mass can be controlled by changing the value of W . 

Finally, we estimate the dissipation coefficients, a and [3. The coefficient a is related to the 
quasiparticle resistance per area rqp as a = l/rqpC'cUp. Here, C = ErEo/d is the capacitance 
per unit area. This relation has been obtained experimentally by Pedersen and Welner.0) 
The quasiparticle resistance rqp obtained below the gap voltage is strongly enhanced at low 
temperatures. Hence, the dissipation coefficient a becomes very small at sufficiently low 
temperatures of mK order. 

To make a comparison, we introduce a dissipation coefficient defined by = l/r-^C'up. 
Here, rn is the normal resistance obtained above the gap voltage. Since is much smaller 
than Tqp, the dissipation coefficient gives the upper limit for a, i.e., a <^ a^- Then, 
can be associated with jc as© 

Jcr„ = ^, (3.10) 

at sufficiently low temperatures {k^T <^ Aq) . Here Aq is the gap of the superconducter at 
zero temperature. From ( p. 3D and ( |3.10|) , we obtain 

an = — (3.11) 
ttAo 



which gives a good estimate for a^- Using Aq = 14[K] and ( p.3|) , we have = 0.088. Thus, 
we evaluate a <^ 0.088. 

The dissipation coefficient (3 also originates from quasiparticle current. It has been observed 
that the ratio between a and (3 is independent of the temperature.^'© Hence, it is expected 
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that (3 is also strongly suppressed at sufficiently low temperatures. We estimate the upper 
bound for (3 as (3 <^ 0.01, which is an experimentally determined value by Davidson et 
a/.Sii* at relatively high temperatures (~ 4[K]). This evaluation for /? is rough, because 
the experimental situation in Ref. |21] is rather different from ones employed in this paper. 
The dissipation amplitude a + /5/3 is, however, not changed drastically even if [3 is enlarged 
several times, because a is expected to be dominant in the dissipation amplitude. 

On the basis of the values given in this section, we estimate the tunneling rate in the later 
sections. 

§4. Tunneling from a Metastable State 

^.1 Tunneling rate 

In this section, we consider the quantum tunneling from a metastable state made by a 
single microresister in a LJJ. First, we formulate the decay rate using the hanger's WKB 
method.lllii) The potential form V{q) is given in (2.45) as 

2e 
cosh^ q 

The ffist term is the driving force due to the external current, and the second term is the 
pinning potential caused by a single microresister. If / is small, the potential V{q) has a 
metastable state at g = go defined by V\qQ) = 0. However, if / is increased and takes a 
critical value fc, the metastable state disappears. The critical value fc is calculated as 

/c = ^ ^ 0.245£(« 1). (4.2) 

When the external current / is taken as / = /c — r/ (r/ > 0), the potential energy has a 
barrier Vq as shown in Fig. 3. To observe quantum tunneling in the laboratory, Vq must be 
small. Hence, we assume 77 ^ /c. The potential is approximated by the quadratic-plus-cubic 
potential around a metastable state as 



V{q) = -2nfq-—^. (4.1) 



where the origin of the coordinate x is located at a metastable state {q = qo + x). From this 
approximated potential, the potential barrier height Vq is calculated as 



Vo = - 12.0s-'/W^'. (4.4) 

The frequency of small oscillation around the metastable minimum defined by uq = 
(l^"(0)/m)^/^ is also obtained as 

c^o=(0)'^'-l.O5.VV^ (4.5) 
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Fig. 3. The potential form V{q) for a fluxon (a) when the external current / is taken as the critical value 
/c, and (b) when / is taken as a slightly smaller value f — fc ^ ri{rj > 0) than fc- 



where the fluxon mass m = 8 has been substituted. Note that Vq and Uq are normahzed 
by Eq and the plasma frequency cUp, respectively. For reference, we consider the exit point 
xo(> 0) defined by V{xo) = 0. It is easily calculated as 



a;„ = W^^ 0.7^/^1/2. (4.6) 

V 87r£: 

The quantum decay rate from the metastable state at T = can be calculated by applying 



Langer's method to the effective action ( |2.8| ). The tunneling rate T takes the form 

T = Aexp{-B). (4.7) 

The exponent B is determined by the action of the nontrivial path qBi^) which minimizes 
the effective action (p.8|). This path is called a bounce solution. As far as Xq <^ 1, Q'bIt) 
satisfies 

kB(r)-gB(r')|«l (4.8) 

for all r, r'. Therefore, the assumption (|2.37| ) is satisfied. The explicit forms of A and 
B are obtained only for the limiting cases. In actual experiments, the damping coefficient 
a + (3/3 in eq. ( p.2|) is small (-C 0.1 in the present situation). Then, the estimate in the weak 
damping limit at T = is obtained by a perturbative treatment a s§) 

/ D X 1/2 

A=VQOiUoiUpi^—j (l + 0(a)), (4.9) 
~ (1 + l-'^4a + 0[^)) (4.10) 
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where a = (a + [3/?>)/2ujq. Since the predominant effect of the dissipation appears in the 
exponent of the decay rate, the correction for the prefactor due to the dissipation is neglected 
in this paper. 

Since takes a small value as given in § ^, the exponent B takes a fairly large value. To 
observe MQT of a fluxon in a time scale of the laboratory, we need to take a value of the 
parameter r] small to increase the decay rate. For example, if we take 77 = 5 x 10"^, e = 0.1, 
then the values of Vq and hujQ are estimated from ( f4.4| ) and ( [4. 51 ) as 

VqEo ~ 0.82W [K] (4.11) 
/icjoc^p ~ 0.34[K], (4.12) 

in the original unit. We note that 77 = 5 x 10~^ and e = 0.1 are not difficult to realize 
experimentally. When W = l[/im] and dissipation effects are neglected (a = /3 = 0), the 
tunneling rate F is estimated from ( [4.7| ) with (4.9) and (4.10) as 

F ~ 2 X 10^[l/s], (4.13) 

which is large enough to observe the MQT in the laboratory. 

4.2 Experimental design to observe MQT of a fluxon 

In this subsection, we propose experimental apparatus to observe the macroscopic tun- 
neling of a single fluxon efficiently. A schematic configuration of the considered Josephson 
junction is shown in Fig. 4. Two superconducting cylinders are separated by a thin insulating 
layer, where one fluxon is captured. The fluxon is accelerated by an externally driven current 
/, which is assumed to be spatially-uniform. A single microresister is made by thickening 
the insulating layer locally. When 1 = 0, the fluxon is captured at the microresister. The 
critical current of the junction is denoted by Jq. 

Experiments on the junctions with a similar topological geometry have already been per- 
formed.i3ii) In these experiments, circular LJJs with the width W ~ 10 [/im] have been 
fabricated. Further, by using scanning electron microscopy, it has been possible to introduce 
individual fluxons into such a system.il' 

The method to observe MQT of a fluxon may be very similar to the experiment technique 
used by Voss and Webb.i* When the external current / = J/Jq is small, the fluxon stays at 
a metastable state caused by the microresister. When the current / approaches the critical 
current fc given by ( [4.2|) , the metastable state vanishes and the fluxon becomes free to move. 
The depinning of a fluxon may occur before / becomes fc because of the quantum tunneling. 
After the depinning, moving fluxons generate a voltage across the junction. Therefore, we 
can probe the depinning of a fluxon by observing an I-V curve of the junction (See Fig. 5) 
Since the tunneling is a stochastic process, the switching current to voltage states has a 



16 



Fig. 4. A schematic configuration of tlie proposed experimental apparatus. An external current / is applied 
to the circular LJJ. Only one fluxon is captured through the junction. A microresister is made in the 
junction. 



distribution. By applying an alternating current, a lot of fluxon- depinning events can be 
observed. (The moving fluxon is captured in repeated intervals which satisfy |/| <^ /c.) The 
probability distribution of the depinning current P{f) is related to the tunneling rate r(/) 
as 

p(/)d/ = r(/) (^1 - J^^ P{f')df'^ dt. (4.14) 

By solving this equation for P{f), we obtain 



Pif) 



df 



dt 



r(/) exp 



df 



dt 



-1 



r(/')d/' 



(4.15) 



where |d//dt| is a sweep rate of the external current, which must be small compared with 
Up. Thus, we may obtain the tunneling rate r(/) by measuring P{f) experimentally.!^ 

The average and the mean-square value of the depinning current are expressed by fc — {rj) 
and Ai] = ijf) — {i])"^, respectively. The values of {i]) and Arj depend on the form of r(/). 
In this section, we only focus on Arj. 

4.3 Simulation of the distribution of depinning current 

To observe MQT of a fluxon in the experiments proposed in the previous subsection, several 
experimental requirements must be satisfled. Among them, requirements for two physical 
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Fig. 5. The expected I-V curve of a circular LJJ. Depinning of a fluxon can be measured by switching to 
voltage states. The averaged switching current and the distribution width are denoted with / — (ry) and 
A77, respectively. The tunneling rate r(/) is obtained by the form of the distribution. 

quantities are important. One is the current accuracy to measure A?7, and the other is the 
temperature T. 

To examine whether the observation of MQT is possible, we first estimate Arj by using the 
estimated values given in § |^. We perform numerical integration of ([4.15|) to obtain P{f)- 
The tunneling rate r(/) is obtained by substituting ( p.8| ) to (4.9) and (4.10) in the original 
unit as 



where Up is estimated in eq. ( |3.3|) , and a = (a + l3/3)/2uo is assumed to be small. Here 
e is taken as 0.1 throughout this paper, which may be realized by thickening the insulator 
layer locally with zero critical current density in a range of O.lAj ~ 2.7 [/im], which can be 
designed in the present available technique. For convenience, we take a serrated wave as 
an alternating current form as shown in Fig. 6. The amplitude of the current is taken as 
twice of fc and the frequency is taken as 1/T = 15 [1/s] to be compared with the experiment 



r(/) ~3.15i/2e^oc^pexp(-5), 

5 ~ 4.1 X 10^r]^/^e-^^^W{l + 1.74a), 
^0 ^ 1.05e'/V^\ 



(4.16) 
(4.17) 
(4.18) 
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Fig. 6. The form of the external current used in the estimate. The amphtude and the frequency of the 
alternating current is taken as 2/c and 1/T = 15 [1/s], respectively. 

performed by Voss and Webb.i^ In this case, we obtain the sweep rate, |d//dt| = 120/c [l/s]- 

Numerical calculation of ( [4 .151 ) for = 1 [/xm] and e = 0.1 is performed for two cases; we 
call the first a strongly dissipative case (a + /5/3 = 0.091), and the second a dissipationless 
case {a = P = 0). Since the dissipation coefficient a + P/3 becomes very small at low 
temperatures in actual experiments, we expect that the distribution of the depinning current 
is almost the same as dissipationless case. To study dissipation effects, however, we also 
consider the strongly dissipative case, in which the damping coefficients are taken as their 
upper bound estimated in § ||. The value of Ar] in the strongly dissipative case gives the 
lower bound for At] in actual experiments. 

The obtained distribution is shown in Fig. 7. From the distribution, we obtain Ar] = 
2.0 X 10"^, {rj) = 2.6 X 10"'* for strongly dissipative case, and Arj = 3.1 x 10"^, {rj) = 4.9x 10~^ 
for dissipationless case. We find that At] is suppressed by dissipation effects. We also find 
that dissipation effects on measurement of the quantum tunneling are not small but modest. 
For the strongly dissipative case, necessary accuracy of the current measurement is given by 
Ar]/ fc ~ 8 X 10"*. On the other hand, in the experiment performed by Voss and Webb,i'' 
AI/Iq < 2 X 10""^ has already been realized. Here AI is the distribution width of the 
switching currents to the voltage states, and Iq is the critical current. Therefore, if circular 
LJJs with W = l[/im] can be fabricated, accuracy of the current measurement to observe 
Arj seems attainable in the present available techniques. 

If we change the junction width W, the distribution width Ar] takes a different value. The 
ly-dependence of Arj is shown in Fig. 8 for both strongly dissipative case and dissipationless 
case. From this dependence, we find that when the junction width W is taken larger, more 
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Fig. 7. The distribution function P{I) with and without damping obtained by numerical integration of 
(|4.15D for jc = 3 X 10^ [A/m^]. The distribution widths are given by A77 = 2.0 x 10"^, (77) = 2.6 x lO""* 
for the strongly dissipative case, and A77 = 3.1 x 10^^, (77) — 4.9 x 10^'' for the dissipationless case, 
respectively. The classical depinning current /c is 0.0245 in this estimate. 



accurate measurement for the current is needed. On the contrary, when W is taken smaller, 
Ai] becomes large, and the observation of MQT of a fluxon becomes easier. In practical 
experiments, we should take a proper value of W so that the measurement of Arj is possible. 
By fitting the estimated data, we obtain Arj oc W~^''^^'^ for dissipationless case. The value of 
the power (~ —0.8) originates from the fact that the exponent is proportional to 'rf/^W. For 
strongly dissipative case, fitting gives Arj oc p^-o-863_ ^qq^i in Fig. 8, dissipation always 
suppresses and the suppression becomes large with the increase of W . This suppression 
is explained by the enhance of the normalized damping strength a = (a + (3 /?>)/2ujq due to 
the decrease of ujq ~ {vV^^- 

Next, we consider finite temperature effects. At high temperatures, thermally- activated 
decay occurs. The decay rate is formulated as 

where Vq is the energy barrier. Here Vq and the temperature T are measured in the original 
unit. It is expected that there exists a crossover temperature T* which separates the thermal 
activated region and the quantum tunneling region. Below T*, the decay rate becomes 
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Fig. 8. The theoretical VF-dependence of Ary with and without damping for jc = 3 x 10^[A/ni^], 
/c = 0.0245. Fitting gives At? oc w~°-^^^ for the strongly dissipative case, and Arj oc W^^-'^^'^ for the 
dissipationless case. 

independent of the temperature. The crossover temperature T* is defined as 



where 7 = a + (3/3, and ciJb is the frequency at the top of the potential barrier. For the 
quadratic-plus-cubic potential, cUb always agrees with uq, which can be controlled by the 
external current f = fc — V- We assume that T* is averaged by the current sweep to the 
value at / = /c — (r/). Then, we obtain the crossover temperature for W = 1 [/xm] as 
T* ~ 54 [mK] for dissipationless cases, and as T* ~ 42 [mK] for strongly dissipative case. 
Here, it appears that the suppression of T* due to damping is not so large. 

To observe MQT of a fiuxon, the junction must be cooled to low enough temperatures below 
T*. We show in Fig. 9 theoretical widths At] expected from the pure thermal activation by 
eq. ( [4.19| ) as well as from the pure MQT with and without damping by eq. (4.16) with (4.17) 
for W = 1 [/xm]. The expected temperature-dependences of Arj are also shown by the bold 
solid curve when both effects coexist. From this calculation, it seems that the crossover 
temperature T* estimated by ( |4.2CI| ) gives a proper criterion. 

Finally, we consider the case that the experimental parameters given in Table I are mod- 
ified. Since the accurate measurement of the parameters, A^, e^, and d is difficult, we only 




(4.20) 
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Fig. 9. The theoretical estimate for temperature-dependence of A77 for the pure thermal activation and the 
pure MQT tunneling with and without damping. The behaviors of A?7 for the case with both effects are 



also shown by the bold thick curves. The crossover temperature T* obtained by ( 4.20 ) is 42 [mK] for the 
strongly dissipative case, and 54 [mK] for the dissipationless case. 



consider jc-control. The critical current density jc can be changed in a wide range by con- 
troUing the thickness of the oxide barrier.!!' The jc-dependence of the estimated quantities 
is given from (|3.2|) , (|373|) and ( |3.5|) by 

Aj cx (4.21) 
Eo^xjI^^ (4.22) 
u^^fj'. (4.23) 

Since g"^ = Tiuj^/Eq is independent of jc, we cannot control quantum fluctuations by jc- 
Instead of it, the damping parameter a can be controlled by jc. When the gap of a super- 
conductor Aq is fixed, jc-dependence of a is obtained from ( |3.3| ) and (|3.11|) as 

a oc j^l (4.24) 

According to this relation, the dissipation effects get important with the increase of jc. The 
jc-dependence of Ar] for W = 1 [/im] with and without damping is shown in Fig. 10. The 
widths At] are almost independent of jc in the dissipationless case, while Arj is suppressed 
in the strongly dissipative case with the increase of jc. 

From (4.23), the crossover temperature T* oc Up is also controlled by jc. The jc-dependence 
of T* calculated from ( [4.20| ) with and without damping for W = 1 [/im] is shown in Fig. 11. 
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Fig. 10. The jc-dependence of A77 with and without damping for W ~ 1 [^m]. 



Fitting the data without damping gives T* oc j^'^^^, which is nearly the value expected from 
(4.23). The deviation originates from the small change in {i]). In the strongly dissipative 
case, T* is suppressed, and the suppression from the dissipationless case grows with the 
increase of jc, because the damping effects become important through ( |4.24| ). At any rate. 



however, the absolute value of the crossover temperature T* is enhanced with the increase 
of jc- Hence, jc should be taken large to observe MQT of a fiuxon more easily in the range 
of jc given in Fig. 11. 

§5. Tunneling in a Two-State System 

In this section, we study a two-state system made by two microresisters. In classical 
mechanics, a fiuxon may stay at either of two stable states at zero temperature. In quantum 
mechanics, however, quantum tunneling through the energy barrier is possible. When the 
energy levels of the ground state at each well is the same and dissipation is neglected, the 
quantum tunneling makes an energy splitting and generates oscillation of a fiuxon between 
the two wells retaining the coherence. This macroscopic effects is called quantum macroscopic 
coherence (MQC). 

In the last ten years, MQC has been investigated theoretically in detaiLi-* MQC is the 
key phenomenon to clarify the validity of quantum mechanics at macroscopic level. 0) MQC 
is far more sensitive to environmental suppression than MQT because a coupling between 
the relevant macroscopic variable and its environment rapidly destroys the phase coherence 
between two states, while we need to retain phase coherence for much longer time for MQC 
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Fig. 11. The jc-dependence of the crossover temperature T* with and without damping for = 1 [/um]. 



than MQT. Thus MQC is in general much harder to observe. Several experiments have been 
proposed to observe MQC on current biased Josephson junctions and SQUIDs. However, to 
date, there is no evidence for MQC in spite of substantial effort. Here, we mainly focus on 
dissipation effects on MQC in LJJs, and pursue possibility to observe MQC in this system. 

When there exists no external current (/ = 0), we obtain the potential made by two 
microresisters as 

cosh^(g-//2) cosh2(g + //2)' 
where the origin of q is taken at the midpoint of the resistors, and the distance between the 
resistors is denoted with /. The potential V{q) has only one stable state for small /, while 
the potential has a double well structure, when I > Iq. The critical length Iq is given by 

/o = In (^^) ^ 1-317. (5.2) 

We only consider the case that the potential barrier is small. In this situation, we may 
assume I — a with a <^ 1. The potential term can be expanded around g = as 

From this potential, the position of the stable states is obtained as q = ±go where 



~ 1.14aV2. (5.4) 
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When a <^ 1, we obtain qq and the expansion in the form of the potential ( |5.3|) is vahd 
around the two stable states. 
The amplitude of the dissipation is determined by a dimensionless quantity K as 



2 



Here m = 8 is the mass of a fluxon. By using the estimated values of g"^, we estimate K as 

~ 3.3 X 10^ + ^ j aVT. (5.6) 

Note that K is independent of details of the potential form and determined by the distance 
between potential minima 2go and the width of the junction W. 

To estimate K, we assume that the ratio between a and (3 is independent of tempera- 
tures.S® Then, because (3 is not dominant in the dissipation amplitude a+l3/3 at relatively 
high temperatures as is observed experiment ally,© we expect that (3 is not dominant also 
at low temperatures. Hence we neglect the /3-term, and estimate the value of K roughly as 

K ~ 3 X lOW— = 3 X lOW— . (5.7) 

Here, we inserted the dissipation coefficient = 0.088 estimated in § The ratio between 
the normal resistance and the quasiparticle resistance r^p can be obtained experimentally. 
The condition to observe MQC is given as 

ii) 

K and ksT < hV/K. (5.8) 

If a = 0.1, ly = 1, and r^/r^p = 10'^ can be realized, the estimated value of i^'(~ 0.03) seems 
to satisfy the first conditions in ( |5.8| ). Although the estimate of the ratio Tn/rqp = 10~^ is 



very rough, it seems to be attainable at low temperatures of mK order.ta^ Hence, we expect 
that the observation of MQC in LJJs appears to be possible within the Caldeira Leggett 
theory. 

§6. Concluding Remarks 

In this paper, we have studied two kinds of quantum tunneling for a single fiuxon. We 
have found that the quantum tunneling from a metastable state may be observed at low 
temperatures (mK order), and that the junction width should be as small as possible. If 
the junction width is taken as ~ required accuracy of current measurement seems 

attainable in the laboratory. We have also found that the observation of MQC in a two-state 
system in LJJs appears to be possible, because dissipation due to the quasiparticle tunneling 
is strongly suppressed at low temperatures. It should be noted, however, that the damping 
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amplitude may increase by other dissipation sources. In the experiment on SQUIDsjHa^ the 
observed damping amphtude is much larger than one estimed by quasiparticle resistance. 
In addition to MQC, we expect other macroscopic quantum effects characteristic in a two- 
state system. For example, it may be possible to observe the incoherent tunneling and the 
population inversion between quantum states in LJJs as observed in SQUIDs.ii'^i) 

The most characteristic feature of LJJs is that the phase difference has a spatial depen- 
dence. In other words, is a field variable with infinite degrees of freedom. This opens the 
possibility of studying combined effects of macroscopic quantum phenomena and many-body 
effects beyond the phenomenological theory by Caldeira and Leggett. In this paper, how- 
ever, dynamics of only one degree of freedom, i.e., the center position of the fiuxon has been 
considered as the first attempt. Other infinite degrees of freedom of the field appears in the 
form of plasmons. They, however, do not play an important role in quantum dynamics of a 
fiuxon because of two reasons: (i) plasmon excitation has an energy gap, and plasmons are 
suppressed at low temperatures; (ii) in the lowest order of g"^, plasmons are decoupled from 
the fiuxon. Hence, within the present study, characteristic features of the field with infinite 
degrees of freedom do not appear. To study many-body effects of the field variable, we must 
consider different situations. Fortunately, it is rather easy to devise LJJs compared with 
other systems described by field theories. New and rich quantum phenomena may appear 
by such devices, and it will be a challenging subject to study the interplay of many-body 
interaction and the macroscopic quantum effects. 
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Appendix A: Derivation of the equation of motion 



In this appendix, we derive the classical equation of motion (|1.2| ) from the model Hamil- 
tonian ( p.l8|) with (|2.11| )-(2.13) within the classical mechanics. First, we write the Hamilton 
equations, 

o(p[x) 

Pj = -T^, (A-2) 

^ 5g (a;) 
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Here 11 = is a momentum conjugate to 0. By substituting the Hamiltonians ( p.ll| )-(2.13) 
into these equations, we obtain series of equations, 

(ptt - + sin0 + / - e6{x) sin0 + —A- + yt?^ = 0' (^-4) 

d(p{x) o(j)[x) 

rrijiij + rrijUj'^qj — Cjcj) = 0, (A-5) 

m'g' + m\ufq\ - c'j(j)^ = 0. (A-6) 



where 

Note that the sign is changed in (A-8) from the integration by part. The Fourier transfor- 
mations of qj, q'j and are denoted by 

<P{x,t) = J^J^kk,u)e^'^^-^-\ (A.9) 
<iAx^t) = j^j^q,{Kuy'^-'-\ (A-IO) 



Fourier transforms of (A-5)-(A-8) and ehmination of qj and q'j give 



5H 

—^ = -iuja{uj)(f){k,uj), (A-12) 



where 



aiuj] 



ik^uj(3{uj)(f){k,uj), (A-13) 



TT 7o cu'(c<;'2 - ^2 - iO+) 
2iu; /■°°, , J/3(cu'l 



/5(cj) = / duj' , , ^ ' -T- (A-15) 

Here, the spectral functions Jaiy^^ and Jpiyj) are defined by (2.14)-(2.15). The positions 
of poles are taken off the real axis to satisfy the causality. When we choose the spectral 
functions as Joiy^^ = Oiu, J/si^) = P^, "we find that a^uj) and P^uj) acquire real parts as 
a'{Lj) = a, 13' {(jj) = j3. Finally the inverse Fourier transforms of (A-12) and (A-13) give 
dissipative terms in the classical perturbed sine-Gordon equation (|1.2|). 
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Appendix B: Details of the Semiclassical Theory 

In this appendix, the semiclassical calculation is presented in detail. Here, we adopt a path 
integral method used by Gervais et a/S The functional integral method obscures problems 
of quantum operator ordering.iil) Actually in the correct formalism based on operator canon- 
ical transformations,^* additional terms, which is absent in the straightforward functional 
approach, appear in the expansion. The extra terms are, however, irrelevant to the present 
calculation because the correction appears only in higher order of g^. 

First we consider the unperturbed sine-Gordon equation given by ( p.26|) and ( p. 27] ). We 
introduce a new field variable (p = (f)/ g to perform the expansion of g^ easily. The partition 
function can be modified as 



Zq = ^ Vip{x,T)VIl{x,T) exp i J dr J dxUip 
JdrJ dx(^^U^ + ^ifl + ^(1 - cosg^)^ 



(B-1) 



where 11 is a field momentum variable conjugate to ip. 
The stationary nontrivial solutions takes the form 



ipo(x — q) = - arctan 
9 



exp(a; — q) 



:b-2) 



where q is the center coordinate of the soliton. Then, we describe the field (p as follows: 

ip{x, t) = ipo{x- q{T)) + r]{x - q{T), r). (B-3) 

Since (p is a. quantum field, it has an infinite number of degrees of freedom in the functional 
integral. Thus, ri{x — q{T),T) also has an infinite number of quantum degrees of freedom, 
and ipolx — q{T)) contains only one quantum degree of freedom in q{T). In order to keep the 
number of degrees of freedom, we should set a subsidiary condition for rj, 



(B-4) 



Here p = x — qi^r) is a soliton-fixed coordinate, and the prime denotes derivative with respect 
to X. In a proper way, we can define momentum variables p(t), 7r(p, r) conjugate to q{T), 
ri{p, r). The partition function is modified as 



Z = ^ Vp{T)Vq{T)Vn{p, r)Vr]{p, t)6{ J dp ^',r])6{ J dp <^[)7r) 
exp ^i J dr pq + i J dr J dpni) — J drH 



:b-5) 



where 



/\2 



H = M, 







(p + / dp 7r?7 
2Mo(l + e/Mo) 



+ 



2 y cosri p 
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- y dp (cosg{(fo + T]) - cos gifo + gr] sin gifo + ]^g'^rf cos g^p^ . (B-6) 
Here Mq = 8/5(^ and 



e = y dp^[,(p)r/'(p,r). (B-7) 
We divide the Hamiltonian into three parts, 

i/.-Mo + ^, (B.8) 

H, = jdp\ (^vr^ + r^'^ + (1 - ^^)^^) , (B-9) 
H^ = H -H^-H2. (B-10) 

Here ifi describes fluxon dynamics in the lowest approximation. The quadratic part H2 
is easily quantized, and gives plasmon excitations. These excitations have a gap of Tiuj^ in 
the original unit. Taking Hi + H2 as the unperturbed part, we can perform the traditional 
perturbation expansion for H^, which includes only higher order terms than quadratic. This 
expansion for (7^ <^ 1 is obtained by Gervais et al^ in detail. The result is 



Zq = Z^j Vq{T)Vp{T) exp i j drp(r)g(r) 
m g'^pi'T'Y 




2m 



+ AE\p{T)] 



(B-ll) 



where m is the soliton mass which equals identically to 8, and Zp is the partition function 
of unperturbed plasmons, which is irrelevant to the tunneling rate. The quantum correc- 
tion AE can be calculated in a systematic way using the traditional perturbation theory 
represented by familiar graphical diagrams. As shown in eq. (B-ll), the lowest order ap- 
proximation for g"^ is nonrelativistic. The Lorentz-invariant form is obtained by summing 
up all relevant diagrams. The result is AE ~ 0{g'^), and there exists no correction to the 
term proportional to p^g^ in AE. 
When g"^ is small, AE can be neglected. After integrating (B-ll) over p, we obtain 

Zo = fVq{T)exp(^-^ Jdr(m+^q{rf^y (B-12) 
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fc - <r|> 



f=l 




logioW [|im] 
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